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SUMMARY 


Research during this period has concentrated on the problem of 
aeroassisted orbital plane change. This maneuver requires the use of three 
impulses - one to deorbit, one to reorbit and one to recircularize at the new 
orbit. The orbit plane change is effected entirely in the atmosphere through 
the use of lift and bank angle control. For circular orbits of nearly equal 
radii, it can be shown that the fuel consumption is minimized by minimizing 
the energy loss in the atmospheric portion of the trajectory. The research 
explores the use of singular perturbation theory to develop an optimal 
guidance law for the atmospheric portion. 

The results to date indicate that singular perturbation methods can be 
applied; however, a difficult terminal boundary analysis is required. The 
reduced solution models only the heading rate dynamics, and produces a 
realistic profile (altitude versus energy) and control to be flown. A large 
terminal boundary layer is required to match the terminal constraint on 
altitude. Most of our effort has been directed at approximate methods for 
solving the terminal boundary layer equations. The equations result from an 
analysis of altitude and flight path angle dynamics on the same time scale. 
A nonlinear control law was derived which produces near optimal results. 
However, the current solution is difficult to implement because it requires 
two switches in the control solution that are heading and altitude dependent. 
In general, the solution is very sensitive to switching times. We propose 
two alternatives to be investigated during the next reporting period. The 
first relies on a linearization of the necessary conditions about the reduced 
solution and the second will examine the analysis of altitude and flight path 
angle dynamics on separate boundary layers. 



1. PROBLEM FORMULATION 


The following three state model has been the subject of our current 
research 


4> = C L pSVAsiny/2mcosY 

(1) 

eh = VsinY 

(2) 

. * 


ey = C^pSV( Acosy + McoSY)/2m 

(3) 


where 

M(h,V) = (2m/C*S)[l-y/V Z r]/pr (4) 

r = r g + h (5) 

and y is the gravitational constant. The objective is to minimize the energy 
1 oss 


J = - E dt 
'o 

where E is the total energy per unit mass 
E = V 2 / 2 - y/r < 0 

The expression for the energy rate in (6) is 
E = -C* (1 + A 2 )pSV 3 /4m 


( 6 ) 


(7) 


( 8 ) 
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where a parabolic drag polar form is used to define the drag coefficient 


C D = C Do + KC L 


( 9 ) 


In the above equations the superscript * denotes the lift and drag 
coefficient values at maximum L/D 


< » <C 0o /K) l/ ’ C* = 2 C Do 


( 10 ) 


The controls are bank angle (p) and the normalized lift coefficient 


x = c L /c L 


( 11 ) 


Note that in this formulation we treat E as constant, but account for the 

energy loss through the performance index. 

In [1] the sensible atmosphere is assumed to occur at h Q = 200,000 ft. 
The starting velocity and flight path angle (V q ,y 0 ) are derived using a 
deorbit impulse aV^ from circular orbit at h c = 100 nm, which is optimized 
for the atmospheric maneuver of interest. The initial heading angle is taken 
as zero. In the SPT formulation, altitude appears as a control variable in 
the reduced problem. The optimal solution has the form 

h * = h(E) (12) 

For comparison purposes, in this study the starting energy is chosen to match 
that of [1], and h Q , V Q are derived from (7) and (12). From conservation of 
energy this results in the same deorbit impulse, but slightly different 
values for h , V . The initial flight path angle is derived from 
conservation of angular momentum. 
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(13) 


Yq = -cos‘ 1 [(r s + h c )(V c - AV 1 )/(r s + h Q )V] 

1/ 

where r g is the mean earth radius and V c = [y/(r s +h)] 2 . The vehicle 

begins the maneuver with a mass m c and, as a result of the deorbit impulse, 
the mass for the atmospheric portion is given by 

m = m c exp (-AV-j/C) (14) 

where C is the characteristic velocity. The terminal conditions are: 

h(t f ) o 200,000 ft, ip(t f ) = * f > 0 (15) 

Since the condition on h(t^) is lost in the reduced solution (12), a terminal 
boundary layer correction is required. 

2. SINGULAR PERTURBATION ANALYSIS 


2.1 Reduced Problem 

Setting e = 0 in (1-3) the necessary conditions for optimality become 


H = A. - E = 0 

0 

(16) 

Y = 0 Acosy = -M 

(17) 

V h o = arg min 

(18) 


h,u 


It can be shown that this results in the following reduced solution: 

A q = (1 + 2M^) V2 (19) 

sinu Q = [(1 + M*)/(l + 2M^)] 1/2 (20) 
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E = const. 


( 21 ) 


h Q = arg 


min 

h 


(v 2 (i 


2 1 / 

M ) /2 } 


where M q is the value of M for h = h Q . The quadrant for the bank angle in 
(20) is resolved based on the following inequalities: 


0 < u Q < tt/2 for M < 0 (22) 

tr/2 < u Q < it for M > 0 (23) 

It can be seen from the above solution that M plays a crucial role in the 
solution process. In [1], M was treated as a constant in the dynamics. 

Since most of the energy is kinetic, V is weakly dependent on h for 
constant E. This can readily be seen from (7) and (5) where changes in h 
give rise to small changes in r. Thus, the minimization in (21) results in a 

value for M very close to zero. The interpretation is that the maneuver 

should be performed at an altitude where gravitational and centripetal forces 
nearly cancel one another. For M small, it can be seen from (19,20) that the 
maneuver is performed at near maximum L/D and at near 90 of bank angle. 
These results are in good agreement with the results in [1]. Figure 1 
compares the altitude profiles derived from (21) with the true optimal 
profile taken from [1], The need for a terminal boundary layer analysis is 
evident in this figure. However, if the vehicle was not required to exit the 
atmosphere, the reduced solution may be sufficiently accurate. 


2.2 Boundary Layer Problem 

A boundary layer analysis is required to obtain a guidance law that will 
both follow the altitude profile defined by (21) (initial boundary layer) and 
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satisfy the terminal constraint on altitude (terminal boundary layer). The 
necessary conditions in the boundary layer are: 


BL 


9H 


where A 


BL 

o 


= + A h VsinT + a y Y " E = 0 

/3L 1 = 0, 9 h BL / 9L 2 = 0 

is determined in the reduced solution from (16) 

= e°/^° 


(24) 

(25) 


(26) 


using the solutions for A Q , y Q and h Q . In (25), and L 2 represent the 
horizontal and vertical components of lift coefficient 


= Asiny L 2 = Acosy (27) 

which are now used as control variables in place of A and y. 

The first condition in (25) results in 
L* = (yVM 1 + M o ) /2 /cos y (28) 


where M , V are the values of M and V corresponding to h = h for the 
o o o 

current value of E. This solution approaches the corresponding reduced 
solution as h approaches h Q . 

The second condition in (25) yields 

L* = -(C*/C*V 2 )A y (29) 

which can also be shown to approach the reduced solution as h approaches h Q , 
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where 


A 0 = C* V ; M /C* (30) 

Y D 0 0 L 

Unfortunately, evaluation of A^ needed in (29) requires the solution of a 
two-point boundary value problem. When close to the reduced solution it may 
be possible to use (30), which results in the following expression for flight 
path angle rate 

^ = C*pSV(Mcosy - \l\ Q /\l 2 )/ 2m (31) 

For y near zero and h near h Q , (31) simplifies to 

Y = C*pSV 0 (M-M q ) (32) 

To obtain a feedback solution for the general case we neglected the 
second term in (24). This was done on the basis that A^= 0 and y is small 
over the entire optimal trajectory. This results in the following explicit 
solution for 

l _2 = - Mcosy + (M cos y - + l)^ 2 (33) 


The first term on the right hand side of (33) is simply the lift required to 
maintain zero flight path angle rate. The second term is always > 0 and 
asymptotically approaches zero as h -*■ h Q and y -► o. Thus this solution also 
asymptotically approaches the reduced solution. Both solutions in (33) 
satisfy the conditions that H gL is minimized and H gL = 0. During the initial 
boundary layer the + sign is used when h < h Q to generate a positive flight 
path angle rate, and the minus sign used when h > h Q . The corresponding 
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value of the costate variable is 


= (Aj * - t )h (34) 

which approaches an indeterminate form (0/0) as h h Q and y 0- The + sign 
is used to initiate the terminal boundary layer. 

At this time repeated trial runs are required to determine the switching 
time so that the desired final heading is achieved when the altitude reaches 
200,000 ft. Also, a characteristic of these profiles is that remains 

close to 1.0 throughout, while M grows to a large negative number near the 
end (on the order of -2.0). This is due to the presence of p in the 
denominator of (4). Thus there is every indication that the sign should be 
switched again in (33) prior to the end of the trajectory so that again 
becomes small. This is also a general characteristic of the optimal profiles 
in [1], 

From (29) it is apparent that should be a continuous function of 
time. There is a discontinuity that occurs at the switch to the terminal 
boundary layer which is a consequence of the singular perturbation 
approximation. A second discontinuity occurs at the second switch which is a 
consequence of neglecting the second term in (24). However, it was observed 
that the second term in (33) passes through a minimum during the ascent 
phase, and the second switch was executed at that time to minimize the 
discontinuity. It is felt that this should more closely approximate the true 
solution if we were able to retain the second term in (24) in the analysis, 
and still preserve an explicit solution for L^. 

A comparison of the resulting flight path with that in [1] for a 40 
plane change is illustrated in Figure 2. Table 1 compares the impulses 
required for the maneuver. Note that the singular perturbation solution 
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Figure 2. Comparison of the guided S.P. solution with the reduced 
solution and the true optimal profile. 
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TABLE 1 


COMPARISON OF TOTAL IMPULSE AND FUEL FRACTION 
REQUIRED FOR A 40° PLANE CHANGE MANEUVER 


GUIDANCE 

LAWS 

DEORBIT 

IMPULSE 

(ft/s) 

BOOST 

IMPULSE 

(ft/s) 

REORBIT 

IMPULSE 

(ft/s) 

TOTAL 

IMPULSE 

(ft/s) 

FUEL 

FRACTION 

OPTIMAL 

125. 

6470. 

177. 

6772. 

.49 

S.P. SOLUTION 

126. 

6642. 

214. 

6982. 

.50 

GUIDED SOLUTION 

374. 

7651. 

122. 

8147. 

.56 

SINGLE IMPULSE 

* 

# 

* 

17497. 

.83 
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results in a fuel fraction close to the true optimal solution, and is 
considerably better than the guided solution in [1]. A comparison to the 
fuel fraction needed for a purely impulsive maneuver is also given which 
clearly demonstrates the advantage of aero-assisted orbital transfer. 
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3. FUTURE WORK 


During the next reporting period we plan to investigate two alternatives 
to constructing a boundary layer solution. The first is based on a 
linearization of the necessary conditions in the boundary layer to obtain a 
linear feedback solution without neglecting the second term in (24) This 
method has been previously used in [2], The second approach analyzes the 
altitude and flight path angle dynamics in separate layers [3]. This 
approach also will yield a feedback solution form, but one which is 
nonl inear. 
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